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Motivated by recent experiments probing anomalous surface states 
of Dirac semimetals (DSMs) NasBi and Cd3As2, we raise the ques¬ 
tion posed in the title. We find that, in marked contrast to Weyl 
semimetals, the gapless surface states of DSMs are not topologically 
protected in general, except on time-reversal-invariant planes of sur¬ 
face Brillouin zone. We first demonstrate this in a minimal 4-band 
model with a pair of Dirac nodes at A; = (0, 0, ±Q), where gapless 
states on the side surfaces are protected only near kz = 0. We then 
validate our conclusions about the absence of a topological invariant 
protecting double Fermi arcs in DSMs using a K-theory analysis for 
space groups of NasBi and Cd3As2. Generically, the arcs deform into 
a Fermi pocket, similar to the surface states of a topological insula¬ 
tor (Tl), and this can merge into the projection of bulk Dirac Fermi 
surfaces as the chemical potential is varied. We make sharp predic¬ 
tions for the doping-dependence of the surface states of a DSM that 
can be tested by ARPES and quantum oscillation experiments. 

Dirac semimetals | Weyl semimetals | Topological insulators | Fermi arcs 

F ollowing the theoretical prediction and experimental dis¬ 
covery of topological insulators mm in the past decade, 
there has been an explosion of interest in understanding the 
role of topology in various quantum states of matter. An im¬ 
portant set of questions concerns the topological properties 
of gapless Fermi systems 0 H 0 El 0 ID- Of particular in¬ 
terest are three-dimensional (3D) semimetals where the bulk 
electronic dispersion exhibits point nodes at the Fermi level 
and the low energy physics is effectively described by Weyl or 
Dirac Hamiltonian 

A striking feature of 3D Weyl semimetals (WSM), which 
necessarily break either time-reversal or inversion symmetry, 
is the existence [5] of topologically protected surface “Fermi 
arcs”. Here the Fermi contour in the surface Brillouin zone 
breaks up into disconnected pieces, which connect the projec¬ 
tion of two Weyl nodes with opposite chirality. The Fermi 
arcs have been recently observed in angle-resolved photoemis¬ 
sion spectroscopy (ARPES) studies of non-centrosymmetric 

TaAs dni El El [13 mils]. 

Here we focus on another outstanding example, the Dirac 
semimetal (DSM), which is the 3D analogue of graphene. In 
the presence of both time-reversal and inversion symmetries, 
the electronic excitations near each node are described by a 
four-component Dirac fermion, and a dispersion that is linear 
in all directions in k-space [T6l [TTI [TS] • ARPES has clearly ob¬ 
served linearly-dispersing bands near Dirac nodes in two DSM 
materials Na 3 Bi [IIEO] and Cd3As2 [21113113]. The signa¬ 
tures of Eermi arcs by ARPES in Na 3 Bi m and by their pe¬ 
culiar quantum oscillations [25] in Cd 3 As 2 [26] have recently 
been reported. Although DSM can in principle appear in a 
system with spin rotational symmetry m, here we focus on 
DSMs in spin-orbit-coupled systems which are more closely 
related to material realizations. 

We can understand the Dirac fermions as two degenerate 
Weyl fermions with opposite chirality, where crystal symme¬ 
tries forbid the two Weyl nodes from hybridizing and opening 
up a gap at each Dirac point [urnsj. Given this picture of the 
bulk, it is natural to expect the surface states in a DSM [TTllTS] 


as two copies of the chiral Eermi arc on WSM surface: i.e. the 
“double Eermi arcs” shown schematically in Eig. [^a). 

In this paper, we address the important question of whether 
the the double Eermi arcs on the surface of Dirac semimetals 
topologically protected? If yes, what is the associated topolog¬ 
ical invariant? The Eermi arcs in WSMs are robustly protected 
for the following reason: each plane that lies between a pair of 
Weyl nodes in momentum space, perpendicular to the separa¬ 
tion between them, has an integer Chern number (5] 1291130] 
associated with quantum Hall effect. The chiral edge modes 
of these momentum-space Chern insulators give rise to robust 
surface Eermi arcs, which end at the projection of bulk Weyl 
nodes. Unlike the WSM, the bulk-boundary correspondence 
in topological phase [8] provides no obvious answer for DSMs, 
since the stability of bulk Dirac nodes require crystal rota¬ 
tion symmetries [171118] which are explicitly broken on open 
side surfaces. It was claimed [25] though, that the surface 
Eermi arcs in DSMs are perturbatively stable against a weak 
symmetry-breaking surface potential, and a strong surface po¬ 
tential can destroy them. 

Here we provide a surprising answer to the question posed 
in the title. We show that the double Eermi arcs on Dirac 
semimetal surface are not topologically protected, and can be 
continuously deformed into a closed Eermi contour without 
any symmetry breaking or bulk phase transition; see Eig. 
The resulting Eermi contours do have topological character, 
even though they are not as exotic as the Eermi arcs in WSMs. 
We show that the surface states of a DSM with an odd (even) 
number of Dirac node pairs are analogous to the surface states 
of a 3D strong TI (weak TI or trivial insulator). 


Significance 

In recent years there has been a surge of interest in quantum 
materials with topologically protected properties that are robust 
against the effects of disorder and other perturbations. An inter¬ 
esting example is the newly discovered three dimensional analog 
of graphene, called Dirac semimetals. These were expected to 
have highly unusual conducting states on their surfaces. We 
show here that the surface states of Dirac semimetals do not 
exhibit the expected topological robustness and, quite generally, 
get deformed into states that are of a rather different character. 
Our theoretical results not only have conceptual importance in 
the field of topological quantum materials but also make clear 
predictions that can be tested in several experiments. 
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To make our results physically transparent, we first focus on 
the 4-band minimal model for DSMs [iTl [181EB and present 
simple arguments for lack of topological protection for surface 
Fermi arcs. We then show numerical results for ARPES spec¬ 
tral functions for the 4-band model which demonstrates how a 
bulk perturbation that respects all symmetries can destroy the 
Fermi arcs. Next, we go beyond the simple low-energy model 
and directly address the two DSM materials of experimental 
interest. We use K-theory [3J [32| [33l [SH [35] to classify the 
topological stability of surface states in the presence of time 
reversal, charge conservation and space group symmetries in 
Dirac semimetals Cd 3 As 2 and NasBi. This classification in¬ 
dicates that topological protection for gapless surface states 
on side surfaces (that do not intersect the £-axis) of Dirac 
semimetals exists only in the high-symmetry plane kz = 0. We 
conclude by examining the experimental implications of our 
results. We discuss how to understand the existing ARPES 
and quantum oscillation measurements in light of our results 
on the lack of topological protection for Fermi arcs, and we 
make sharp predictions for testing our conclusions. 

Minimal model 

We begin with a simple 4-band model for a Dirac semimetal 
with 4-fold rotational symmetry along the z-axis defined by 
the Hamiltonian: 

H(k) = e^-{-^{coskx-\-cosky — 2)-\-tz{coskz — cos Q)]r. 

-\-XsmkxCrxTxXsinkyayTx. [1] 

The Pauli matrices a act on spin and r in orbital space, t and 
tz are hopping amplitudes and A is spin-orbit coupling (SOC). 
Ignoring for the moment, there is a gap everywhere in the 
Brillouin zone (BZ) except at two points k = (0, 0, ±Q) where 
the low energy spectrum is described by linearly-dispersing 
Dirac fermions. 

We can always add to eq. M an that preserves all 
the symmetries and vanishes at me Dirac nodes, e.g., = 

ti{cos kz — cos Q)-\-t 2 {cos kxcosky— 2) or = ti{coskz — 
cos Q) {cos kx-\-cosky). While this term does not qualitatively 
change the Dirac spectrum in the bulk, it gives rise to a cur¬ 
vature of the Fermi arc in the surface BZ, as discussed below. 

The symmetries of H (k) and their representations in terms 
of spin and orbital operators are as follows. Time-reversal 
is implemented by the antiunit ary operator 0 = icr^ • /C, 
where JC is complex conjugation and 0^ = —1. Inversion I 
corresponds to Ui = r^, with the two orbitals having even 
and odd parity. Two-fold rotation C 2 ^x{y) about the x (or y) 
axis is implemented by Uc 2 ^{y) ~ and n-fold rotation 

Cn,z by = exp (i7rcr2(l — 2T2)/n), where n = 4 here. 

if(k) also has mirror reflection symmetries with respect to 
the j — a:, y, z-planes, with kj —kj, which is implemented 
by Ur. = Uc 2 jUi = icTjTz. 

Instability of surface Fermi arcs 

Surfaces with a normal not parallel to z-axis are expected to 
support gapless states with Fermi arcs, shown schematically 
in Fig. [Ra) for side surfaces. We address the lack of topolog¬ 
ical protection for the Fermi arcs of H (k) from three different 
points of view, (i) First, we examine the stability of ID (one¬ 
dimensional) edge states that exist on 2D slices of BZ at a 
fixed, generic value of kz- Here we define a generic kz to be 
kz ^ {±Q,0,7r}, so that a fixed-Zc^ plane will neither contain 
Dirac nodes, nor will it have the additional symmetry arising 
from time-reversal invariance, (ii) Next, we prove the exis¬ 


tence of an extra mass term in the gapped Dirac Hamiltonian 
that describes the plane with a generic kz- (iii) Finally we 
show how adding fully symmetric, bulk perturbations to the 
3D Hamiltonian H{\c) can destroy the Fermi arcs, without de¬ 
stroying the bulk Dirac nodes. The projection of these bulk 
nodes onto the surface does remain gapless, as do the 2D slices 
located at = 0 as they have extra (time-reversal) symmetry. 

Without loss of generality we consider a (100) surface. This 
surface preserves only the following symmetries: time re¬ 
versal 0, mirror reflections Ry,Rz and their combinations. 
Therefore a fixed-plane in a system with with a (100) 
surface only respects two symmetry operations: mirror re¬ 
flection URy = icFyTz w.r.t. y-plane, and the combination 
0i?2 = OUr^ = icTxTz • JC of time reversal and mirror w.r.t. 
i-plane. Note that Gr^ = 1. Under these symmetries the 
Hamiltonian at fixed kz, denoted by H{kx,ky), transforms as 

G^lH{kx,ky)GR^ = H{-kx,-ky) and H{kx, ky)URy = 

H{kx, J^y) • 

(i) First consider the surface BZ {ky,kz) for a (100) sur¬ 
face. If we look at a ID slice at fixed-Zc^ where we might 
have expected gapless edge states described by Hedge = 

- kF)ipky Here = {tpR,ky,i’L,ky) is 
a two-component spinor of right (R) and left (L) movers and 
jl are Pauli matrices in (R, L)-space. Note that the dispersion 
vf and the location of the Fermi arc kp are both, in general, 
dependent. 

The two symmetry operations on the edge states are repre¬ 
sented by Gr^ — ijjLx ' JC and URy — i/iy. We now add per¬ 
turbation M = '^kyl^y'^ky, wMch preserves both these 

surface symmetries: [GR^,yy] = [URy,yy] = 0. Clearly M is 
a mass term for Redge that destroys the gapless edge states. 

(ii) Next, we gain insight into the instability of surface states 
on a fixed-plane within classification scheme [32l [33] [34] of 
topological insulators. For a fixed, generic kz, our system is a 
gapped 2D insulator described by a Dirac Hamiltonian 

Hd = kx^x + ky'jy + myo, [2] 

up to an additive constant term. Here ^x = ctxTx, 7 ^ = cTyTx 
and 7 o = Tz are Dirac matrices in the minimal model. Let 
us ask if we can add another mass term m '70 to Hd which 
preserves all the symmetries and anti-commutes with all the 
Dirac matrices 7cc,^,o- If we can do this, the m > 0 phase can 
be continuously deformed into the m < 0 phase without clos¬ 
ing the gap at m = 0 , and Hd would be topologically trivial. 
Such a term indeed exists: 70 = cfzTx commutes with Gr^ 
and URy and anticommutes with 'yx,y,o- Therefore the surface 
states at generic kz are not topologically protected. 

In contrast to generic kz, the time-reversal-invariant planes 
kz = 0 and kz = tt are gapped 2D insulators with the well- 
known Z 2 topological index associated with quantum spin Hall 
(QSH) effect In the DSM materials Cd 3 As 2 and Na 3 Bi, 

there is a nonsymmorphic glide reflection gc containing half¬ 
translation along i-axis, which serves as a mirror reflection 
for kz = 0 , 7 r planes satisfying {gc)‘^ = —1(+1) for kz = 0('7r). 
Generically such nonsymmorphic glide reflections guarantees 
that the Z 2 index must be trivial for the kz = tt plane [36] . 
which will be elaborated below. Therefore in these materials, 
only kz = 0 plane can support a nontrivial Z 2 index. A sim¬ 
ple calculation Ea shows that the kz = 0 plane is indeed a 
nontrivial QSH insulator, and gapless surface states at = 0 
are topologically protected. 

(iii) Finally we directly examine the gapless surface states 
of the 3D Hamiltonian 0 in a slab geometry with a sur¬ 
face at X = 0 , which is translationally invariant along y and 
z. We numerically calculate the spectral density A(k,ct;) = 
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— (l/ 7 r)ImG(k, ct; + iO^) associated with the electron Green 
function G(k, uj). Here we measure uj with respect to the zero 
chemical potential, which coincides with the bulk Dirac nodes. 
We show in Fig. [^the surface [panels (b,c,d)] contributions 
to H(k, 0) as a function of in-plane momentum (ky,kz)- Note 
that ARPES probes precisely this spectral density multiplied 
by the Fermi-Dirac function. 

In Fig. Wb), we plot the result for H(k.) of eq. 0 and 
clearly see the two Fermi arcs extending from one Dirac node 
to the other. In the absence of the term of the form dis¬ 
cussed above, H(k.) has an (unnatural) “chiral symmetry” 
and the two Fermi arcs become “degenerate” and collapse to 
ky = 0 for —Q <kz<Q [38]. With in place, we get the arcs 
shown in Fig. I^b). 

Now we ask n there is a perturbation that has the full sym¬ 
metry, does not shift the bulk nodes and yet has a (Ja G 
structure that anticommutes with each term in H(k.) (other 
than the term, of course). Such a perturbation does exist 
and it must be of the form 

( 5 i/ 4 (k) — m!(coskx — cosky) sinkz CzTx [ 3 ] 

This perturbation destroys the Fermi arcs near the projections 
of the Dirac nodes at kz = on to the surface BZ as shown 
in Fig. We clearly see this in panel (c) where ( 5 i/ 4 (k) with 
m' = OAt deforms the Fermi arcs into a closed Fermi pocket. 
An increase in the strength of the perturbation to m' = 0.8t 
shrinks the arcs to even smaller pockets as shown in Fig. md). 
While the surface Fermi arcs are progressively destroye(i7the 
bulk Dirac nodes (and bulk states near Ep = 0) are unaf¬ 
fected by these perturbations. We note that this perturbation 
of 0(k^) is higher order than that retained in usual k • p per¬ 
turbation theory. 

The minimal 4-band model m is equivalent, up to a unitary 
transformation [38| , to the well-known effective k • p Hamilto¬ 
nian around the T point for DSM materials nziiiH]. Despite 
its simplicity, this model with n = 4-fold rotational symmetry 
about the z-axis captures the band inversion near the T point 
of BZ in Gd 3 As 2 . A similar model can also be adopted for the 
case of NasBi with n — 6 -fold rotational symmetry [38j . 


K-theory 

Using K-theory we next show that the surface states and asso¬ 
ciated Fermi arcs in Gd 3 As 2 and Na 3 Bi are not topologically 
protected. This mathematical framework has the great ad¬ 
vantage of focusing only on the space group, U{1) charge con¬ 
servation and time reversal symmetry of a material, without 
the need for a low-energy effective model such as the one used 
in the analysis above. We thus obtain general and rigorous 
results that are directly applicable to real materials. 

In the K-theory approach [32l (33] [34] , the classification of 
distinct gapped symmetric phases is reduced to the follow¬ 
ing mathematical problem: what is the “classifying space” 
of symmetry-allowed mass matrix for a generic Dirac Hamil¬ 
tonian preserving certain symmetries? Different symmetric 
phases correspond to disconnected pieces of the classifying 
space, which cannot be continuously connected to each other 
without closing the bulk energy gap. Mathematically the 
group structure formed by these different phases is given by 
the zeroth homotopy group 7 ro( 5 ') of classifying space S. Here 
we only sketch the idea behind the calculations; detailed anal¬ 
ysis are given in Supplementary Materials [38] . 

Na 3 Bi: The hexagonal DSM material Na 3 Bi has a unit 
cell is defined by a = a(l,0,0), b = a(—1/2, \/3/2,0) and 
c = c(0,0,1). Its space group is centrosymmetric P 63 /mmc 
(No. 194) [H], with its symmetry operations are generated 


by 6 -ford screw sq\ (x,y,z) {x — y,x,z 1 / 2 ), glide re¬ 
flection Qc: {x,y,z) (y, 1 / 2 ), inversion /: {x,y,z) 

(—X, —y, —z), and Bravais lattice translations Ta^h.c- 

The two Dirac nodes HD in Na 3 Bi are located at = 
±(5o(27r/c)(0, 0,1). We consider two different surfaces which 
can potentially contain nontrivial Fermi arcs, (i) A (010) sur¬ 
face spanned by the a and c axes, and (ii) a ( 110 ) surface 
spanned by (a — b) and c. In each case we proceed as follows. 
We begin by determining the reduced set of symmetries that 
are preserved in the presence of the surface. By analyzing how 
these symmetries act on a general Bloch Hamiltonian, we de¬ 
termine those symmetries that preserve a 2D plane in the BZ 
zone at a fixed value of kz ^ All such planes, except 

those passing through the nodes at are fully gapped. 

The basic strategy is to exploit the classification of gapped 
insulators described by Dirac Hamiltonians defined on these 
2D planes embedded in the 3D Brillouin zone. We use the 
K-theory classification Ha ESI E] of gapped free-fermion 
quantum phases by examining the combined algebra for the 
Dirac 7 -matrices and the symmetry operations and determin¬ 
ing whether the classifying space for symmetry-allowed mass 
terms is connected or not. 

Using this approach we conclude that a generic 2D plane in 
the BZ at fixed kz is topologically trivial and does not need to 
support protected surface states. The plane kz = tt/c is also 
shown to be trivial, due to the presence of the nonsymmoprhic 
glide reflection Qc. Only the kz — ^ plane is found to have a 
nontrivial Z 2 QSH index. 

Cd 3 As 2 : The space group of the tetragonal DSM Gd 3 As 2 is 
centrosymmetric I4i/acd (No. 142) [39]. In addition to Bra¬ 
vais lattice translations, its symmetry operations are gener¬ 
ated by 4-fold screw rotation along z-axis: {x,y, z) (y, 1/2 — 
a:,l/4 + z), glide reflection gc’. {x,y,z) (x, —y, 1/2 + z), 

and inversion I: {x,y,z) {—x,l/2 — y, 1/4 — z), where 

( 54 )^ = (gc)^ = —Tz. The Bravais lattice is expanded by 
T±i,±i,±i = (±l,±l,±l)/2 similar to BGG lattice. The two 
Dirac nodes in Gd 3 As 2 are located along the i-axis. 

The K-theory procedure is used to determine the topolog¬ 
ical nature of various 2D planes in the BZ at fixed kz, in a 
manner completely analogous to that described above, taking 
into account of the symmetries of Gd 3 As 2 . The calculations 
are described in supplementary materials [38] , and the conclu¬ 
sions are identical to the ones given above. Only the kz = 0 
plane has a nontrivial Z 2 QSH index, all other fixed kz planes 
are trivial. 


Nature of Surface States 

We next show that the closed Fermi pocket on a side surface 
is essentially the same as the surface Dirac fermion of a 3D 
strong topological insulator (TI). Gonsider a weak perturba¬ 
tion that breaks Cn symmetry, but preserves time reversal 
and inversion. We require the strength of this perturbation is 
smaller than the bulk gap at the time reversal invariant mo¬ 
menta (TRIM) of the 3D BZ. Glearly this perturbation will 
gap out the Dirac nodes but will not affect of parity eigenvalues 
of filled bands at the TRIM. Now that the bulk is fully gapped 
band insulator with time reversal symmetry, its Z 2 topolog¬ 
ical invariant is given by the number of odd-parity Kramers 
pairs in filled bands at TRIM as shown by Ref. m Since 
the kz = 0 plane is a 2D QSH insulator and kz = tt is n 
trivial 2D insulator, the gapped bulk must be in a 3D strong 
TI phase. In the process of turning on the perturbation, the 
surface states at = 0 remain gapless since time reversal is 
always preserved. This demonstrates the equivalence between 
the closed Fermi pocket for DSMs and surface Dirac fermions 
of strong TIs. 
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An important conclusion from our analysis is that if the 
bulk band touching occurs at an odd number of pairs of Dirac 
nodes not located at TRIM (such as in Cd 3 As 2 and NaaBi), 
then the surface states cannot be fully removed due to time 
reversal symmetry. Let Np Dirac points be located in the top 
half of the bulk BZ (0 < kz < tt) and their Np time-reversal 
counterparts in the bottom half BZ (—tt < kz < 0). Based on 
our arguments, this will lead to Np closed Fermi pockets on 
a generic surface, which are either centered around = 0 or 
kz = TT. If there are an odd number of Fermi pockets centered 
around = 0, the gapped kz — ^ plane in 1st BZ must cor¬ 
respond to a 2D QSH insulator with an odd number of helical 
edge states. Therefore, if Np is odd, kz — ^ and kz — tt must 
have opposite Z 2 topological indices, with one plane being a 
trivial 2D insulator and the other a QSH insulator. Hence the 
surface states of of an odd-Np DSM cannot be fully gapped 
out, in complete analogy with a 3D strong TI. On the other 
hand, if Np is even, both the kz = 0 and kz — tt planes share 
the same Z 2 index i.e. they are both trivial 2D insulators or 
both QSH insulators. In the former case there are generally 
no surface states, while gapless surface states in the latter case 
can be gapped out by translation symmetry breaking pertur¬ 
bations in analogy to a 3D weak TI. As a result, an “odd 
DSM” with an odd number of pairs of Dirac points must sup¬ 
port time-reversal-protected gapless surface states, in contrast 
to an “even DSM”, whose surface states can either be fully 
gapped or need protection from crystal translation symmetry. 

Experimental Implications 

We have shown above that, in general, the surface states of 
Dirac semimetals Cd 3 As 2 and Na 3 Bi have a closed Fermi sur¬ 
face (FIG. (c)), in contrast to the double Fermi arcs that 
one might naively expect (FIG. (b)). We now address the 
important question of how one can experimentally distinguish 
between these two types of surface states. 

Assume, first, for simplicity that the chemical potential is 
at the Dirac point; (we discuss below the case when it is not). 
The qualitative differences between the two outcomes are then 
as follows. The double Fermi arc has a singular change in slope 
when two arcs meet at the projection of a Dirac node on to 
the surface BZ. In contrast, a Fermi contour has a smooth 
curvature everywhere in the surface BZ and it does not pass 
through the two points that are the projections of the bulk 
Dirac nodes. In addition, the wave-function associated with 
the surface states merges with the bulk at the tip of the arc, 
namely the Dirac node; in contrast the states associated with 
a regular Fermi contour are exponentially localized at the sur¬ 
face. All of these features could, in principle, be used to 
distinguish between the two scenarios using sufficiently high- 
resolution ARPES data with photon energy dependence used 
to separate bulk and surface contributions. 

We next consider the very interesting quantum oscilla¬ 
tion experiments on Gd 3 As 2 thin films m- The observed 
Shubnikov-de Haas oscillations seem consistent with “Weyl 
orbits”, which consist of the combination of bulk Landau lev¬ 
els (LLs) and surface Fermi arcs [25]. How can we understand 
these experimental results, if the surface states have a closed 
Fermi pocket instead of double Fermi arcs? 

First of all, both transport gomilES] and ARPES [42l[22] 
experiments confirmed that Cd 3 As 2 is slightly n-doped, in 
the sense that Eermi energy Ep lies above the Dirac points 
with kp — 0.04A . Therefore the surface projection of the 
bulk electronic states, near each Dirac point, is in fact a small 
“blob” denoted by a red oval in EIG. ^ (d)-(f). In this case, 
the surface states at the Eermi energy may (EIG.[^f)) or may 
not (EIG.[^e)) be connected with the bulk Eermi blobs, de¬ 


pending on Ep. When the surface pocket merges into the bulk 
blobs, the magnetic orbits responsible for quantum oscillations 
must be the Weyl orbits proposed in Ref. |25j which consists 
of both bulk LLs and surface Eermi arcs. Shubnikov-de Haas 
oscillations in Ref. [2^ suggest that this is the situation in 
Gd 3 As 2 . The available ARPES results in Na 3 Bi M are also 
consistent with surface states merging into bulk projection as 
kz increases. 

However if we start to introduce p-type dopants to the sys¬ 
tem and lower the Eermi level, ultimately the surface pocket 
must get disconnected from the bulk blobs, as shown in EIG. 

(d)-(e), when the Eermi level is close enough to the Dirac 
points. In this case the surface Eermi pocket will provide a 
closed 2D magnetic orbit for cyclotron motion. This pocket 
is analogous to the 2D Dirac fermion on the surface of topo¬ 
logical insulators (TIs) and will lead to quantum oscillations 
similar to those observed in TIs [43j |44j [45], with the elec¬ 
tron acquiring a Berry phase of tt as it goes around the Eermi 
contour. As with any 2D Eermi surface, the quantum oscil¬ 
lation frequency Es has a 1/cosO dependence on the angle 0 
between magnetic field and surface normal direction. There is 
an important difference between the quantum oscillations in 
the two scenarios where the bulk and surface states are mixed 
and separated is that in the latter case, Es has no dependence 
on the thickness of the sample, in sharp contrast to the Weyl 
orbits [25] [26] in the former. 

Therefore quantum oscillation experiments at different dop¬ 
ing levels provide a sharp distinction between double Eermi 
arcs and Tl-like closed surface Eermi pockets, serving as an 
experimental test of our conclusion. There will be two fre¬ 
quencies in quantum oscillations on DSM thin films [251126] : 
Eh associated with bulk Eermi surfaces (due to deviation of 
Ep from Dirac point) and Es related to the surface states. 
As we move Ep towards bulk Dirac points by doping the sys¬ 
tem, Eh will always decrease monotonically. As demonstrated 
in Ref. [26], a triangle-shaped sample does not exhibit quan¬ 
tum oscillations with frequency Es , while a rectangular sample 
does, because the Weyl orbits depend on the thickness of the 
system along the field direction. Thus in a triangle-shaped 
sample, double Eermi arcs cannot lead to quantum oscillation 
at Es independent of the Eermi level. On the other hand, as 
the surface pocket is disconnected from the bulk blobs by tun¬ 
ing Ep close to the Dirac point (see EIG. [He)), the frequency 
Es, with a 2D angular-dependence, will show up even in a 
triangular sample. In particular, when the magnetic length 
Ib — \/h/eB satisfies kd ^ where kd is the shortest dis¬ 
tance between the surface pocket and the bulk blob in surface 
BZ, the closed surface pocket cannot tunnel into bulk LLs 
to form a closed Weyl orbit. When the Eermi level lies ex¬ 
actly at the Dirac point, kd ~ 0.04A“^ gives an estimate of 
B < lO^T which is within the experimental reach. Of course, 
very high-resolution ARPES studies can also directly reveal 
the separation of surface Eermi pocket from the bulk blobs as 
the photon energy is varied. 

Concluding remarks 

Let us conclude by summarizing our main results. Each bulk 
Dirac node looks like two copies of a Weyl node, which naively 
suggests that the surface electronic structure of a DSM should 
look like two copies of that of the Weyl semimetal (WSM), 
namely a double Eermi arc. We show here - using a simple 
4-band model and rigorous K-theory calculations - that the 
double Eermi arcs are not topologically protected, unlike the 
Eermi arc in a WSM. We find that an arbitrarily small bulk 
perturbation, which preserves the Dirac nodes and all the sym¬ 
metries, can lead to a deformation of the double Eermi arc into 
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a more conventional Fermi contour. Nevertheless, we show 
that the surface states in a DSM cannot be completely de¬ 
stroyed, because there must be topologically protected states 
on the kz = 0 plane, similar to a 3D TL 

Finally, we explore in detail the experimental implications 
of our results. We show how we can reconcile our results on 
the lack of topological protection with existing experiments 
on NasBi [24] and Cd 3 As 2 [26]. While these have been inter¬ 
preted as being consistent with Fermi arcs, we argue that they 
are actually a consequence of the mixing of bulk and surface 
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Fig. 1. (Color online) (a) Schematic k-space picture of a Dirac semimetal showing Dirac nodes along kz axis in bulk Brillouin zone (BZ) and possible double Fermi arcs 
on the surface BZ’s, shown as blue squares. Note that surfaces perpendicular to 2 : axis have no arcs. A 2D slice of the bulk BZ perpendicular to kz axis is shown as a green 
square, which projects to a green dashed line on a side surface, (b) Surface spectral density of model in which clearly shows the existence of double Fermi arcs on (100) 
surface. (c,d) Continuous deformation of double Fermi arcs on the (100) surface by adding the perturbation (51^4 (k) to (m. (c) shows the effect at m' = —OAt, while (d) 
corresponds to m' = —0.8t, showing that the Fermi arcs are progressively destroyed by increasing strength of perturbation. The red solid circles in (b-d) correspond to the 
projection of bulk nodes and we set Ep = 0 to line up the Fermi level with bulk Dirac nodes. (e,f) correspond to electron doped systems with m' = —0.8t (case in d) 
by raising the Fermi energy to Ep = O.lt and Ep = 0.15t, respectively. The large red blobs in (e,f) mark the projection of bulk states onto the surface BZ. The Fermi 
contour of the surface states is disconnect from the blobs in panel (e), while it merges into the bulk states in panel (f). 
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